
A SINGLE FIELD OF VIEW METHOD 
FOR RETRIEVING TROPOSPHERIC 
TEMPERATURE PROFILES FROM 
CLOUD-CONTAMINATED RADIANCE DATA 


Donald B. Hodges 


Prepared by 

TEXAS A&M UNIVERSITY 
College Station, Tex. 77840 
jor George C. Marshall Space Flight Center 




NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


AUGUST 1976 


tech library kafb, nm 



TECH LIBRARY KAFB, NM 


1. REPORT NO. 2. GOVERNMENT ACCESSION NO. 

NASA CR-2726 


& TITLE AND SUBTITLE 

A Single Field of View Method for Retrieving Tropospheric 
Temperature Profiles from Cloud-Contaminated Radiance Data 


ODblMSb 


3. RECIPIENT'S CATALOG NO. 


5. REPORT DATE 

August 1976 


6. PERFORMING ORGANIZATION CODE 


7. AUTHOR(S) 

Donald B. Hodges 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 

Center for Applied Geosciences 
Texas A &M University 
College Station, Texas 77840 

1 2 SPONSORING AGENCY NAME AND ADDRESS 

National Aeronautics and Space Administration 
Washington, D. C. 20546 


8 . PERFORMING ORGANIZATION REPORT # 

M-174 


10. WORK UNIT NO. 


11. CONTRACT OR GRANT NO. 

NAS8-26751 


13. TYPE OF REPORT & PERIOD COVERED 


Contractor 


14. SPONSORING AGENCY CODE 


IS. SUPPLEMENTARY NOTES 


Prepared under the technical monitorship of the Aerospace Environment Division, Space 
Sciences Laboratory, Marshall Space Flight Center, Alabama 


16. ABSTRACT 


A method is presented to retrieve single field of view tFOV) tropospheric temperature profiles directly from 
cloud-contaminated radiance data through the use of auxiliary data such as observed shelter temperature and 
estimated cloud-top height. The iterative technique utilized is an extension of the work of Chahine as modified by 
Smith and Duncan. It Is shown that a well-defined temperature profile may be calculated from the radiative transfer 
equation (RTE) for a partly cloudy atmosphere when the average fractional cloud amount and cloud-top height for the 
FOV are known. A cloud model Is formulated to calculate the fractional cloud amount from an estimated cloud-top 
height. 

The method is then examined through use of simulated radiance data calculated through vertical Integration of 
the RTE for a partly cloudy atmosphere using known values of cloud-top height(s) and fractional cloud amount(s). 
Temperature profiles are retrieved from the simulated data assuming various errors in the cloud parameters. 

Temperature profiles are retrieved from NOAA-4 satellite-measured radiance data obtained over an area 
dominated by an active cold front and with considerable cloud cover. Excellent radiosonde data for 11 stations 
participating in the Atmospheric Variability Experiment (AVE ill) available for the area are used for comparison 
with temperature profiles retrieved from the NOAA-4 data. The effects of using various guessed profiles and the 
number of Iterations are considered. 

The results of the Investigation indicate: 

(1) The single FOV method presented Improves the accuracy of guessed profiles when retrievals from 
cloud-contaminated radiance data are attempted. It 1 b significant that in the method presented. Improvement in 

the guessed temperature profile was noted under the cloud layer where retrievals using other single FOV techniques 
tend to deteriorate. 

(2) Through proper choice of a guessed profile, It appears possible to Improve many profiles Independent 
of the cloud amount present. 

(3) Use of the method presented should provide useful data for mesometeorologlcal research In the form 

of an accurate temperature profile representative of a relatively small horizontal area, or many temperature profiles 
to show the pattern of temperature change over a larger area. However, the absolute accuracy of the retrieved 
profile Is also a function of the a priori knowledge of the state of the atmosphere possessed by the researcher. 

(4) It Is possible to make an accurate estimate of the average tops of a thick overcast layer through use of 
the cloud model when there are no thick clouds present above the overcast layer. 

(5) Observed cloud parameters do not possess sufficient precision for use directly in the RTE to retrieve 
temperature profiles. 


18. DISTRIBUTION STATEMENT 



Category 47 


19. SECURITY CL ASSIF. (of thU r. party 20. SECURITY CLAS SI F. (of this p««.) 21. NO. OF PAGES 22. PRICE 


Unclassified 


Unclassified 


$5.25 


For sale by the National Technical Information Service, Springfield, Virginia 22161 
















TABLE OF CONTENTS 


Page 

Abstract i 

TABLE OF CONTENTS v 

LIST OF FIGURES vii 

LIST OF TABLES xi 

LIST OF SYMBOLS AND ACRONYMS xii 

1. INTRODUCTION 1 

a. Statement of the problem 1 

b. Objectives 3 

2. THEORETICAL DEVELOPMENT 4 

a. The Radiative Transfer Equation for a_ clear 

atmosphere 4 

b. Related cloud models 12 

3 . A NEW METHOD FOR RETRIEVING TEMPERATURE PROFILES FROM 

THE RADIATIVE TRANSFER EQUATION FOR A PARTLY CLOUDY 
ATMOSPHERE 22 

a. General 22 

b. The method 22 

c. The cloud model 25 

4. PARAMETRIC STUDY 29 

a. Procedures 29 

b. Results 30 

5. APPLICATION TO AVE III DATA 45 

a. Procedures 45 


iii 



TABLE OF CONTENTS (Continued) 


Page 


b. Data 45 

c. Results 51 

d. Applicability to mesometeorological 

research 81 

e. Possible sources of error 82 

6. CONCLUSIONS , 83 

7. RECOMMENDATIONS FOR FURTHER RESEARCH 86 

REFERENCES 88 


lv 



LIST OF FIGURES 


Figure Page 

1(a) VTPR Transmittance Curves (McMillin et al. 1973) . . 6 

1(b) VTPR Weighting Functions (McMillin et al. 1973) . . 7 

2 Temperature sounding retrieved from Duncan's 

Method for a clear atmosphere using simulated 
radiance measurements for broken (0.80) low cloud 
conditions 31 

3 Temperature sounding retrieved from Duncan's 
Method for a clear atmosphere using simulated 
radiance measurements for scattered (0.18) high 

cloud conditions 32 

4 Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere. Simulated radiance 
measurements were prepared for broken (0.80) low 

cloud conditions 33 

5 Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere. Simulated radiance 
measruements were prepared for scattered (0.18) 

high cloud conditions 34 

6 Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere accomplished by em- 
ploying a guessed profile exhibiting sharp 

temperature inversions 35 

7 Temperature sounding retrieved from the RTE for 

a partly cloudy atmosphere accomplished by em- 
ploying in the retrieval the same values of the 
cloud parameters as used to calculate the mea- 
sured (simulated) radiance values ......... 37 

8 Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere accomplished by em- 
ploying in the retrieval the same cloud-top heights 
as used to calculate the measured (simulated) rad- 
iance values but fractional cloud amounts that 

are each 0.1 less than the values used to calcu- 
late the simulated radiances 38 

9 Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere accomplished by em- 
ploying in the retrieval the highest cloud-top 
height used to calculate simulated radiance mea- 
surements and calculating a one level fractional 

cloud amount (0.27) at that level . 39 


v 


LIST OF FIGURES (Continued) 


Figure 

10 

11 

12 

13 

14 

15 

16 
17 


Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere accomplished by em- 
ploying in the retrieval the highest cloud-top 
height used to calculate simulated radiance mea- 
surements and calculating a one level fractional 
cloud amount (0.59) at that level 

Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere accomplished by em- 
ploying in the retrieval a cloud height that is 
significantly lower than the highest cloud-top 
height used to calculate simulated radiance mea- 
surements and calculating a one level fractional 
cloud amount at the significantly low level . . . . 

Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere accomplished by em- 
ploying in the retrieval a cloud-top height that 
is significantly higher than the highest cloud- 
top height used to calculate simulated radiance 
measurements and calculating a one level frac- 
tional cloud amount at the significantly high 
level 

Temperature sounding retrieved from the RTE for 
a partly cloudy atmosphere. It was accomplished 
by calculating a value of fractional cloud amount 
at an estimated cloud-top height, employing the 
calculated value to retrieve a temperature pro- 
file and then calculating revised values of cloud 
amount and temperature until calculated radiance 
values and simulated measurements converge 

Satellite radiance measurements (o) at 0233 GMT 
and AVE III radiosonde runs (•) at 0000Z, 6 
February 1975 

Rawinsonde stations participating in the AVE 

III experiment (Fuelberg and Turner, 1975) 

Surface synoptic chart for 0000 GMT, 6 February 
1975 (Fuelberg and Turner, 1975) 

Centerville, Ala. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 
III radiosonde data 


Page 

40 

41 

42 

44 

47 

48 

49 

59 


vi 


LIST OP FIGURES (Continued) 


Figure Page 

18 Jackson, Miss. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 60 

19 Shreveport, La. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 61 

20 Stephenville, Tx. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 62 

21 Del Rio, Tx. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 63 

22 Midland, Tx. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 64 

23 Nashville, Tenn. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 65 

24 Little Rock, Ark. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 66 

25 Monette, Mo. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVE 

III radiosonde data 67 

26 Amarillo, Tx. 6 Feb 1975 retrieved temperature 
profile compared with the guessed profile and AVe 

III radiosonde data 6*3 

27 Marshall Space Flight Center, Ala. 6 Feb 1975 
retrieved temperature profile compared with the 

guessed profile and AVE III radiosonde data .... 69 

28 AVE III 839 mb isotherms and surface frontal 

positions 70 

29 Retrieved 839 mb isotherms and surface frontal 

positions 71 

30 AVE III 699 mb isotherms and surface frontal 

positions 72 

vii 


LIST OF FIGURES (Continued) 


Figure 

31 

Retrieved 699 mb isotherms 
positions 

and 

surface 

« • • « 

frontal 


Page 

73 

32 

AVE III 509 mb isotherms and surface frontal 
positions 

• • 

74 

33 

Retrieved 509 mb isotherms 
positions . 

and 

surface 

frontal 

» • • » • 


75 

34 

AVE III 412 mb isotherms 

» • 

» • • • « 

» • • • • 


76 

35 

Retrieved 412 mb isotherms 

• ■ 

• • • • 

* m • • • 


77 

36 

AVE III 299 mb isotherms 

. . . 

. . . . 

» • • • • 


78 

37 

Retrieved 299 mb isotherms 


> • • • 4 

i • • • • 


79 


vm 


LIST OF TABLES 


Table Page 

1 Estimation of average cloud-top heights 50 

2 RMS errors (°C) below cloud-top level for 

various guessed profile . 54 

3 RMS errors (°C) for retrieved profiles 56 

4 RMS errors for guessed (G) and retrieved 

profiles 57 


IX 



LIST OF SYMBOLS AND ACRONYMS 


U 

A* 


Total 


B 

B 

r 

CO, 

c 

/\ 

c 


c , c 
12 


H 2° 


I* 


used to denote the power to which weighting function 
raised 

effective cloud amount (A=Ne) 

effective cloud amount for the lower of two layers 

effective cloud amount for the upper of two layers 

a function of A T and A 
L u 

sum of effective cloud layers 
Planck function (radiance) 

Planck radiance at frequency r 
carbon dioxide 
a numerical correction 
a computed estimate of C 
constants used in computing temperature 

the quantitative difference between the average radiances 
arising from the clear and cloudy portions of a field of 
view (FOV) 

water vapor 

clear column radiance 

measured radiance for a cloudless atmosphere 
measured cloud-contaminated radiance 


I* a computed estimate of I* 

I measured radiance 

MEAS 

I average radiance arising from the cloudy portion of the FOV 

^CLR average radiance arising from the clear portion of the FOV 

(equivlanet clear column radiance when obtained by filtering 
cloud effects) 

I clear column radiance averaged over many FOVs 


x 



r 


x 


max 

I CLR' 


if j 
k 


1 


m 

n 

N 

N* 


O 


3 


P 


P 

o 


P 

c 


u 


r 

T 

T 

o 

T 

c 

W 


03 

03 


LIST OF SYMBOLS AND ACRONYMS (continued) 

a calculated estimate of clear column radiance 

i m ^ n maximum and minimum possible values of clear coluinn 

CLR , . 

radiance 

radiance arising from a black-body cloud 
subscripts 

the number of effective cloud layers 

the number of FOV combinations used to calculate average 
clear column radiance 

the number of frequencies employed to calculate temperature 
superscript used to denote iteration number 
fractional cloud amount 

ratio of fractional cloud cover in two FOVs 

ozone 

pressure 

surface pressure 

pressure at cloud-top level 

pressure at the top of the lower of two cloud layers 
pressure at the top of the upper of two cloud layers 
residual; error 

a nominal frequency near the center of the CO^ band 

temperature 

surface temperature 

cloud-top temperature 

weighting function 

window channel frequency 

a short-wave window channel frequency 


xi 


LIST OF SYMBOLS AND ACRONYMS (continued) 

o) a long-wave window channel frequency 

L 

X, Y functions of cloud- contaminated radiance used to obtain 
equivalent clear column radiances from measured values 

2/7 

x a function of pressure (x = P ) 

x value of x at the surface 

o 

x^ value of x at cloud top level 

P 

c 

Ab change in the Planck function 

Ab weighted average of changes in the Planck function 

e emissivity 

e cloud emissivity 

c 

ri a function of the fractional cloud amount in the two fields 

of view 

y mass of absorbing gas above a given level 

v frequency 

v average frequency in a small interval 

t transmissivity 

tco^ transmissivity of carbon dioxide 

tH 2 ° transmissivity of water vapor 

"to^ transmissivity of ozone 

AGL above ground level 

AVE Atmospheric Variability Experiment 

FOV field of view 

HIRS High Resolution Infrared Sounder 

ITPR Infrared Temperature Radiometer 

JPS Joint Planning Staff 

MSFC Marshall Space Flight Center 


xii 



LIST OF SYMBOLS AND ACRONYMS (continued) 


NASA 

NOAA 

RMS 

RTE 

VTPR 


National Aeronautical and Space Administration 

National Oceanic and Atmospheric Administration 

root mean square 

radiative transfer equation 

Vertical Temperature Profile Radiometer 


xiii 



I. 


INTRODUCTION 


a. Statement of the problem 

Over the past decade, a number of studies have dealt with vari- 
ous aspects of the problem of retrieving temperature profiles from 
satellite-measured radiances in the infrared region of the electro- 
magnetic spectrum. Much of the research in this area has had the 
ultimate purpose of providing profiles that are suitable for use as 
input to numerical forecast models . 

Comparatively little attention has been given to the possibil- 
ity of using satellite-derived profiles for mesometeorological 
research. Originally, this was probably because of the poor spatial 
resolution of the radiance measurements of the early satellites. As 
satellite technology has advanced, the spatial resolution for a 
single FOV has improved. However, the most accurate temperature 
profiles retrieved from cloud -contaminated radiance data have gen- 
erally been achieved through use of a multiple FOV method in which 
cloud filtered and clear radiance data obtained for many FOVs are 
averaged . 

Through use of a multiple field of view cloud model (the term 
"cloud model" will be used to represent techniques of filtering or 
otherwise accounting for clouds in the RTE) an average profile 
generally representative of a large horizontal area is retrieved. 

For the High Resolution Infrared Sounder (HIRS) of the Nimbus 6 
satellite a resolution of 200 km is obtained (Smith eit al_. 1975) . 

Recently, good results also have been obtained through use of 
measurements in the microwave region which are relatively unaffected 
by cloud contamination. For the Nimbus 5 satellite, average kinetic 
temperatures were obtained for layers of approximately 10 km thick- 
ness centered near 4, 11, and 18 km (Staelin, 1974) . Temperatures 
at discrete levels are then obtained through use of the correlation 



between these temperatures and the temperature over the weighting 
function layer (Waters et al . 1975) * The Nimbus 6 satellite con- 
tains a microwave sensor with a resolution of 145 to 330 km from 
nadir to scan limit that will have a maximum cloud-caused error of 
2K over water and IK over land (Staelin et^ al . 1975) . The imple- 
mentation of microwave techniques and improvements of multiple FOV 
cloud models have not altered the fact that the cloud problem re- 
mains the most serious obstacle in the retrieval of temperature 
profiles for mesometeoro logical research. 

In the present investigation an attempt is made to provide 
useful retrievals for this purpose by improving the spatial re- 
solution of the temperature profile through use of a single FOV 
cloud model based on observed cloud and temperature data. The 
cloud and temperature data are used in conjunction with real and 
simulated radiance data from the NOAA series of satellites that are 
currently used for operational retrievals over ocean areas where 
cloud amounts are not too great to retrieve significant information . 
Previous single FOV models have been based on climatology (Rodgers, 
1970) or, as in the Smith et_ al. (1970) model, have been found most 
useful in reducing the influence of clouds on the solution profile 
above cloud-top level (Fritz et^ al . 1972) . 

Routinely observed cloud amounts and heights have not been 
used as input to the cloud models used in temperature retrieval 
work. The unknown emissivity of the observed clouds and the un- 
certainty of cloud-top heights, fractional amounts, and number of 
cloud layers are some of the reasons why these data have not been 
utilized. Fritz et al. (1972) pointed out some of the difficulties 
involved in using auxiliary cloud data to assist in determining 
the effective cloud cover (Ne) , height, and amount, but conceded 
"...this has not been tried yet, so it is not known what effect 
such a procedure would have on the accuracy of temperature retriev- 
als." For most models in current use the cloud-contaminated data 
are filtered to obtain an "equivalent clear column radiance" prior 
to solving the radiative transfer equation (RTE) to retrieve the 
temperature profile. 
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This research represents an extension of previous research in 
that a method is devised to use auxiliary cloud information to ob- 
tain an improved single FOV temperature profile directly from the 
RTE for a cloudy atmosphere. Also, in order to implement the pro- 
posed method a cloud model is devised to obtain the fractional cloud 
amount at an estimated cloud-top height through use of a search 
among radiance values which are calculated at a single frequency 
for various cloud amounts . A single FOV method provides much 
better spatial resolution than a multiple FOV method, and the method 
presented in this research has yielded improved temperature profiles 
below cloud layers. This is the region in which retrievals from 
other single FOV techniques tend to deteriorate. 

b. Objectives 

The objectives of this research were to: 

(1) Develop a theoretical model to retrieve single FOV 
temperature profiles from cloud-contaminated radiance data? 

(2) Examine the model through use of a parametric study uti- 
lizing simulated radiance data computed from a known temperature 
profile to investigate the errors in retrieved profiles caused by 
errors in cloud amounts, cloud heights, and other parameters; and 

(3) Compare the profiles retrieved through use of the model 
against radiosonde profiles obtained during the Atmospheric Var- 
iability Experiment (AVE III) . 

The procedures used to achieve the above objectives will be 
fully explained in the sections that follow. 
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2 . 


THEORETICAL DEVELOPMENT 


a. Ihe Radiative Transfer Equation for a_ clear atmosphere 

Over the past decade attempts to determine atmospheric temper- 
ature profiles have concentrated on the 15-ym carbon dioxide band 
of the electromagnetic spectrum. To retrieve temperature profiles 
the RTE may be solved numerically for frequencies in this band. 

Later in this section methods for solving the RTE will be discussed. 
However, the following fundamental assumptions are required in all 
methods of determining temperature profiles from carbon dioxide 
band data (JPS, 1973): 


(1) The mixing ratio of CO^ is constant below 30 km. 

(2) The atmosphere below 50 km is in local thermodynamic 
equilibrium (i.e. Planck's function and Kirchoff's Law may be used). 

(3) Scattering by aerosols is negligible. 

As a practical matter it is also necessary to assume that 
v = v (defined below) in the small spectral interval of each channel. 

Using the above assumptions, the RTE for a plane parallel, 
cloudless, and non-scattering atmosphere may be expressed as 


I (v) = B[v,T(P 0 ) ] • t[v,y (P Q ) 3 - 



B [v , T (P) ] 


(1) 


where I(v) is the radiance (intensity), B[v,T(P)] is the Planck 

radiance, and / t[v,y(P)] is the transmittance at frequency v of the 

mass of absorbing gas y above pressure P, and the partial derivative 
9 t[v,y (P) 1 

^ is a weighting function giving the relative atmospheric 

contributions to I(v). A set of integral equations which is used 
to calculate radiances at the various frequencies of the 15— ym 
band may thus be obtained from a known temperature profile when 
T.[v,y(P)] and thus can b e calculated. 

There are several methods of calculating values of t[v,y(P)3 
used to evaluate the weighting functions. In this research the com- 
puter program used in the calculations is the same as that used for 
the NOAA series of satellites. A general discussion of the procedure 
used for each gas with references to original sources is given by 
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McMillin et al. (1973). Transmittance s for carbon dioxide (t__ ), ozone 

CC>2 

(t_ ) , and water vapor (t _) were calculated by use of the following 

H2O 

equation 


tfV/y(P) 1 


= t C o 2 tv ^ (p)] 


[v,y(P)]-T [v,y(P)l 


(2) 


to obtain the total transmittance above a given pressure level. The 
ozone transmittance is a relatively minor correction to the total 
transmittance at a given level. However, as the distribution of 
moisture is highly variable it is necessary to use a guessed moisture 
profile as input to the computer program. The mixing ratio of carbon 
dioxide is relatively constant with height, and thus estimates of the 
variation of transmittance with pressure have been determined. 

Transmittance and weighting function curves for the 15-ym chan- 
nels of a Vertical Temperature Profile Radiometer (VTPR) instrument 
(McMillin et al. 1973) are shown in Fig. 1. 

For each channel (labeled 1-6 in Fig. la) the values of trans- 
mittance above a given atmospheric level increase as lower values of 
pressure are used to indicate the atmospheric level until they 
approach 1.0 asymptotically at the top of the atmosphere. This 
result should be anticipated. As the mass of CC> 2 present at an 
atmospheric level decreases with an increase in height above Earth's 
surface, the transmittance of radiation emitted at a given level 
must increase with height above the surface. 

As transmittance increases with increasing height above the sur- 
face, the change of transmittance with height (i.e. the weighting 
function) increases to a maximum value and then decreases until the 
value 0.0 is approached at the top of the atmosphere (Fig. lb) . 

Thus for a given channel the contribution to the measured radiance 
according to (1) from a given level will increase with altitude 
until the peak on the weighting function is reached and then de- 
creases above. The height at which the weighting function peaks 
for a particular frequency is dependent on the location of the 
frequency with respect to the center of the 15-ym band. At fre- 
quencies close to the centers of absorbing bands a small amount 
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C0 2 TRANSMITTANCES 



Fig. 1(a). VTPR Transmittance Curves 
(McMillin et_ al ♦ 1973) . 
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PRESSURE (mb) 


C0 2 WEIGHTING FUNCTIONS 



Fig. 1(b) . VTPR Weighting Functions 
(McMillin et al. 1973) . 
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of the absorbing gas will attenuate much of the transmitted radia- 
tion and, therefore, most of the outgoing radiation near the centers 
of absorbing bands arises from the upper levels of the atihosphere 
(JPS, 1973) . At frequencies far from the centers of absorbing bands 
it takes a large amount of gas to attenuate much of the transmitted 
radiation; thus, most of the outgoing radiation away 1 from the 
centers of absorbing bands originated in the lower levels of the 
atmosphere . 

Due to an overlap in the weighting function curves (Fig. lb) 
the amount of independent information about temperature that can be 
obtained from radiance measurements is limited. After a certain num- 
ber, providing additional radiance measurements for more channels 
in the 15-pm band will be redundant in the sense that the integral 
equations will no longer be independent. Seven degrees of freedom 
is considered a limit. Also, because of the shape of the weighting 
function curves, small variations in measured data caused by instru- 
ment error may lead to large errors in the final solutions of the set 
of equations and, therefore, to an unrealistic temperature distribu- 
tion (JPS, 1973) . 

Measured radiances are dependent on air temperature, instrument 
characteristics,, and atmospheric transmittance. The distribution of 
CC^ in the atmosphere is assumed constant and this does not introduce 
serious errors. However, water-vapor content and the distribution 
of some gases other than CO^ is quite variable and atmospheric 
transmittance can be significantly affected (Fritz et al . 1972) . 

Early studies by Kaplan (1959) and Yamamoto (1961) revealed the 
necessity of a special stabilizing technique in the computation of 
the temperature profile. Wark and Fleming (1966) suggested a practi- 
cal method of overcoming the instability problem, based on the expan- 
sion of the deviation of temperature from standard or mean values in 
terms of orthogonal functions. A regression method based on the 
statistical relationship of the distribution of temperature and radi- 
ance measurements was successfully derived for operational use by 
Smith et al. (1970) . 



Iterative methods have also been developed. Iterative methods 
work well when there is a reasonable guessed profile (Allison et al. 
1975) . Significant changes in the retrieved profile compared with 
the guessed profile identify a poor guess . The "minimum information" 
iterative method (Smith et al. 1972) is currently used to retrieve 
temperature profiles from the Infrared Temperature Profile Radiometer 
(ITPR) of each Nimbus series satellite, and until March 13, 1975 
(Werbowetzki , 1975) was used to retrieve temperature profiles from 
the Vertical Temperature Profile Radiometer (VTPR) of each NOAA 
series satellite. 

An iterative method of retrieving temperature profiles was form- 
ulated by Chahine (1968, 1970) . Upon using the mean value theorem 
and the fact that the Planck function has a stronger dependence on 
temperature than the weighting function, the relaxation formula 

B[v. ,T (n+1) (P) ] = Blv. ,T (n) (P) ] (3) 

1 l Kn) (V ± ) 1 

is derived from the RTE. In the above formulation, I (v . ) is measured 
( n ) 

radiance and T ^ is the temperature obtained on the nth iteration. 

A temperature profile is retrieved through use of (3) as follows. 
First, an initial guessed profile of temperatures corresponding to 
the approximate peaks of weighting functions is devised. Then, em- 
ploying the guessed profile, numerical integration of the RTE is 
accomplished and a radiance value is calculated for each of the 
sounding frequencies . If the residuals 



I(v) 


approach zero for the. individual frequencies then the guessed tempera- 
ture profile is a solution. When convergence is not obtained, a new 
guess for the temperature values corresponding to each of the i 
sounding frequencies (one temperature per frequency), T^ n+1 ^(P^), is 
required. A combination of the relaxation equation above with the 
mathematical expression for Planck's Law for the given frequencies 
leads to 
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(5) 


T i°(pJ) = c 2 v i /ln '{ 1 -f 1 ~ ex P(c 2 Vi/ T i (n) (P±) )] (I i <n) /S i )j 


Numerical integration is again accomplished with subsequent iteration 
until convergence is obtained. Chatline’ s method is non-linear and 
should converge for a wide spectrum of guessed profiles. However, as 
pointed out by Barcilon (1975), Chatline’ s method is only valid for 
a square matrix of frequencies and levels. Consequently, only a 
limited number of solution points (for temperature) may be retrieved 
through use of this method. 

Smith (1970) subtracts the iterative form of the RTE 


i (n) (v.) = B[v.,T (n) (P q )3 t(v.,P 0 ) + 
( B[v.,T (n) (P) 1— 9 ( - p i/P - dP 


( 6 ) 


from (1) and through use of the assumption that B[v^,T^ n+1 ^ (P) ] - 
B [v^, T ^ (P) ] is independent of pressure over the sensed atmospheric 
layer obtains the iterative equation 


B[v.,T (n+1) (P) ] = B(v, ,T (n) (P)] + [I (v . ) 
1 


x (n) (V.)] 


(7) 


Here, the radiance values -are computed by numerical integration from 
guessed temperature values that are not restricted to the number of 
the sounding frequencies employed. When convergence of measured 
and calculated radiances is not obtained the equation 


< n+1 > (V..P) = C-V./ln c 1 v 3 + BW.,T (n+1 >(P) 1 


2 r 


B[v i ,T (n) (P)] 


(8) 


is then used to calculate independent estimates of the entire temp- 
erature profile for each frequency. To obtain the best estimate of 
temperature at each level, T^ n+1 ^(P), from the independent tempera- 
ture estimates for each frequency, the weighted average 
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T (n+1) 


(P) 


E 


T (n+1) 


in 

(V ± ,P)W(V i# P)/ 

i=l 


W(V i ,P) 


where W(v^,P) = d t{v^,P) if P ^ P Q 
and WCv^Pq) = t(v i ,P 0 ) 


(9) 


(i.e. the W values for each frequency are the weighting functions 
used to compute radiance from the RTE) is computed for the m frequen- 
cies . 

Duncan (1974a) observes that substitution of (3) into (6) is 
equivalent to multiplying the iterative form of the RTE by I(v^)/ 
i' ’ (vj and hence the computed Planck function satisfies the RTE 
exactly . 

He then uses Chahine's relaxation formula to implement Smith's 

concept of obtaining an independent temperature profile estimate from 

2/7 

each radiance measurement. Upon scaling pressure values by x = P 
[following Smith et_ al_. (1972) ] for pressure levels from 1000 mb to 
0.01 mb (for accuracy in applying the trapezoidal rule to the RTE) 
and substituting (3) into (8) , independent estimates of temperature 
are given by 


(n+1) 


3 t (n) 


, , x , CiV. I (V.)/I(V.) + B[V.,T 

(v^x) = c v ± /ln 1 1 1 1 1 


(n) 


(x)] 


( 10 ) 


Brv i ,T (n) ( x )] 


for each frequency at each of 100 pressure levels. A weighted average 
temperature similar to that obtained from (9) may then be computed. 

A comparison of various aspects of the minimum information. 

Smith, and Duncan methods, is given by Alexander (1974.) 

Chow (1975) demonstrates that use of the weights of (9) , which 
are the weights of the integral form of the RTE, makes it impossible 
to determine any of the fine structure of the atmosphere that is not 
present in the initial guessed profile. Furthermore, the retrieved 
profile must retain a shape similar to the guessed profile unless 
the weights are raised to some power, a. He concludes that 
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increasing the value of a not only increases vertical resolution, but 
also increases the rate of convergence. However, these results are 
achieved at the expense of a significant increase in the effects of 
radiance measurement errors on the retrieved profile (i.e. computa- 
tional instability) . Chow, therefore, concludes that a small value 
of a should be used in determining profiles for Earth's atmosphere 
since climatological and forecast profiles are available for use 
as guessed profiles. 

b. Related cloud models 

The presence of clouds causes serious complication and results 
in errors in retrieved temperature profiles . In order to demonstrate 
this fact, simulated measured radiance data (radiances calculated 
through vertical integration of the RTE from a known temperature 
profile) were calculated for various cloud amounts and tops . Re- 
trievals were then attempted from this data using Duncan's method. 

No attempt was made to modify the procedures used for a clear atmo- 
sphere through use of a cloud model. Examples of results will be 
shown later. Attempts to account for the effects of clouds have 
results in models proposed by Smith et al. (1970), Rodgers (1970), 
Chahine (1970), and Jastrow and Halem (1973), among others. The 
Rodgers ' model is based on the correlation between temperatures 
above and below clouds . The others are based on the equation 


I MEAS NI CD 


+ (1 - N)I 


CLR 


( 11 ) 


where I is the measured radiance, and ^ are the average 

MEAS CD CLR 

radiances arising from the cloudy and clear portions of the FOV, and 
N is the amount of cloud cover. Through use of radiance measure- 
ments in a cloud-dependent channel, two or more channels and/or 
FOVs, (11) may be transformed to a set of simultaneous equations and 

solved for N. The known N is then used in (11) to obtain I which 

CLR 

is then used as the measured radiance in the temperature retrieval. 
With improved resolution of radiance elements a multiple FOV tech- 
nique of eliminating cloud cover has proven successful in improving 
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the accuracy of retrieved profiles for use in numerical weather pre- 
diction. In general, however, satellite-derived temperature profiles 
have not provided a significant input to mesometeorological research 
and the accuracy of profiles derived from any model yet formulated 
tends to decrease with increasing cloud cover. 

Below are outlined some of the important models that have been 
developed to filter cloud-contamination effects from the measured 
radiance data. 

(1) Smith (1968) Model . It is assumed that average radiance 
arising from two resolution elements (radiance spot measurements) is 
the same; thus implying that cloud heights are the same if the clouds 
are black-bodies. The FOV of the sensing system must be as small as 
practicable to insure no radical difference in cloud heights between 
elements. It is also assumed that angular resolution is sufficiently 
high that each of the resolution elements encloses an area much 
smaller than the area for which the average temperature is desired. 

Measured radiance for frequency is given by (11) . For two 
different elements, subscripts 1 and 2, (11) may be rewritten 


and 


W^V - N i I CD 1 (V i ) + (1 - V^LR/V 


( 12 ) 


■’•MEAS^V ^2 J CD 2 + (1 N 2 )I CLR 2 (V i ) * 


(13) 


If ^ N 2 , these equations can be solved for average clear column 

radiance. If I„„,„ ^ for the window channel (in this chan- 

MEAS^ MEAS 2 

nel transmittance is nearly 1.0 for atmospheric gases, but not for 
cloud), then cloud cover for the two elements is not the same. For 
small and adjacent elements it is assumed that 


I„ (v . ) = l__ (v. ) 
CD ' l CD, 1 

2 1 


(14) 


and 


I CLR (V i ) I CLR 1 * V i ) I CLR 2 (V i ) 


(15) 
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With the above assumptions in mind, (12) and (13) are solved for 
I CLR (V i ) t0 glVe 


I CLR^ V i' ^MEAS, (V i' " N * I MEAS. 


(v t ) 3/(1 - N*) 


( 16 ) 


where 


N* = N x /N 2 


(17) 


When N is not equal to N , and I 


MEAS- 


is not equal to I 


MEAS2 


(i.e. 


only differences in fractional cloud cover cause variations in mea- 
sured radiance) , N* can be obtained from two simultaneous measure- 
ments of radiance in the window (m) region through use of (12) and 
(13) . Clear-column radiance may then be calculated through use of 
(16) . If instead it was desired to calculate N* from clear column 
radiance obtained from a surface temperature ^ observation or measured 
directly by a high angular resolution window radiometer, the following 
equation was used: 


N* = [I 


MEAS, 


^ I CLR^'* // ^MEAS. 


(w) - 


I CLR^ ) ] 


(18) 


In any case, if the field of view of the sensing element is small, 
average clear column radiance may be calculated over a large geo- 
graphic area from the equation 

WV = I VCLR. j/ L <19) 

j=*l j=l 

where W^ are the weights (1-N_.*) . These weights are used because 

the observations tend to be inflated by — rr [Ref. (16)] . For the 

( 1— N * ) 

342 different combinations of adjacent elements in a 10 x 10 matrix 
of spatially independent observations originally used with this 
model (i.e. £=342), the effects of random observational errors and 
differences in cloud heights were assumed to be reduced to 
insignificance. 

(2) Smith et al. (1970) Model . A first guess of the temperature 
profile is made in order to compute corrections to the observed radi- 
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ances. Then using radiance measurements in two spectral intervals 
sensitive to clouds, a system of two equations in two unknowns is 
formulated and solved for the equivalent clear column radiance. The 
equivalent clear column radiance is then used to make a next guess 
temperature profile. The iteration is continued until convergence is 
obtained. In this method the temperatures above clouds may be 
improved, but temperatures below are highly dependent on the first 
guess (Fritz et al . 1972) . 

This is a two level model. It is based on the "...common atmo- 
spheric situation..." where a senri transparent layer of cirrus exists 
over opaque middle clouds. It can be shown that the following equa- 
tion describes the effect of cloud cover for an atmosphere containing 
no more than two layers of cloud: 


W v) - 


I MEAS ^ 


A X [V,P ,P , T (P) ] 

U. U L 


+ A*Y[v,P l ,T(P) ] 


(20) 


where I (v) is the radiance that would be measured in channel (v) 
LLR 

under clear conditions ("the equivalent clear column radiance"), 

I..™., (v) is the measured radiance in channel (v) , and A , called the 
MKAb U 

fractional radiative cloud amount (effective cloud cover) for the 
upper cloud layer, is the product of the fractional cloud amount and 
the fractional cloud transmittance for the upper layer. Also 


A* = A t (1-A ) (21) 

L u 

where A is the fractional radiative cloud amount for the lower cloud 
Li 

layer. Further, the parameters X and Y are given by 


and 


X = B[V,T(P ) ] t(V,P ) - B[v,T(P )]t(V,P ) 
s s u u 



p 

u 


B[V,T(P)]dt(V,P) 


(22a) 
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(22 b) 


Y = B[V,T(P )]t(V,P ) - B[V,T(P t ) ] t(V,P_) 
S S L L 



B[v,T(P)3dt(v,p) 


/ 


respectively. The essence of (20) is that a correction 


C(V) = A^X + A*Y (23) 

must be computed and added to the measured radiance to get the equiv- 
alent clear column radiance. In computing C(v), clouds are allowed 
to exist at any two "standard" pressure levels below 150 mb. An 
estimate of equivalent clear column radiance from the radiative 
transfer equation is first computed from a guessed temperature profile. 
This was done for each of the channels most sensitive to clouds for 
Nimbus 3. An estimate of C(v) is therefore 


I CLR^ V ^ ~ J MEAS^ 


(24) 


for v=714, 750, and 899 cm 1 (window), and where ^ CLR ^) is th e esti- 
mated clear column radiance calculated from a guessed temperature pro- 
file. Next, the X and Y terms of (22) are specified from the guessed 
profile for all standard pressure-level combinations. Substituting 

C(v) of (24) for C (v) in (23), A and A* are calculated for all pos- 

u 

-1 -1 

sible standard pressure-level combinations for 714 cm and 899 cm 
channels by solving the simultaneous equations generated by the sub- 
stitutions. Using the various pressure-level combinations, different 

values of A and A* (and thereby A) are used to calculate C(v) from 
U -1 L 

(23) for the 750 cm channel. The value of C(v) is then calculated 

-1 

from (24) using 750 cm . The most probable cloud cover is specified 
as that for which 


C (v) - C (v) | = Min 


(25) 


where v is 750 cm 


As new temperature profiles are calculated, the 
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cloud cover computations are repeated until the cloud correction 
ceases to change from one iteration to the next. It should be 
noted that in this model the surface temperature is specified in the 
guessed profiles to prevent them from becoming unrealistically cold 
(corrections are always from cold to warm) . Also, this method is 
generally useful only in determining temperature profiles down to 
cloud-top level, and profiles below clouds were predominantly cal- 
culated from statistical relationships between the temperatures above 
the cloud layer and those below. A variation of the above method has 
been recently formulated by Cooper (1975) 1 , but as yet satisfactory 
results have not been obtained. 

(3) Rodgers (1976) . In this method Rodgers uses the high cor- 
relation with temperatures above cloud level to obtain temperatures 
below. He also suggests that the accuracy of this method can be 
improved by using other data sources such as surface temperature, 
cloud picture data, and forecast profiles. The basic approach is 

to obtain a maximum probability estimator of the atmospheric state 
(Fritz et. al . 1972) . 

(4) Smith et al. (1974) Multiple FOV Approach . This is the 
method used by NOAA to filter ITPR (Nimbus 5) data. Two geographi- 
cally independent observations are required. But the observations 
must be close as it is assumed that temperature profiles and, 
therefore, equivalent clear column radiances are the same for both 
observations . Error is introduced if measured radiance variations 
are caused by anything but variation in cloud amounts. A good esti- 
mate of surface temperature is also required in this method as is a 
high resolution measuring device (Fritz et al. 1972) . When variation 
of radiance from one FOV to another is due to variation in cloud 
amount only and surface temperature is known, then clear column 
radiances may be computed from two sets of independent data. 

It is first necessary to determine if cloud properties (height, 
opacity, etc.) in the FOVs chosen are similar. It should be noted 
that the Nimbus 5 ITPR has window channels at both 3.7 ym and 11 ym. 


Cooper, M., 1975: Personal communication. 
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The same brightness temperature would be observed for both channels 
for a uniform and opaque scene (e.g. Earth's surface), but brightness 
temperatures would differ significantly for the two channels under 
broken cloud conditions because of the different dependence of the 
two channels on temperature. Where clouds are present their pro- 
perties are considered the same for the two FOVs if 


min 

CLR 


( “s> 


20 £ < 

1-N* 


I CLR,i, j <“s> 


< 


max 

CLR 


(co ) + 
s 


2oe 

1-N* 


(26) 


where 


min 

I 

CLR 


and I 


max 

CLR 


are the minimum and maximum possible clear column 


radiance values [Ref. Smith et al . (1974) for the necessary proce- 
dures to determine and w refers to the short wave (3.7 ym) 

CLR CLR S 

window channel and the term added to (subtracted from) I (I ) is 

CLR CLR 

the expected error of I . . (go ) , the clear column radiance for the 

CLRf X f J S 

two FOVs. Also, 


N* = N i = ri (co ) - I . .(go)] / 

N i,i — L MEAS,i L CLR, l,] L 

N. 

3 

^MEAS, j ^L^ I CLR,i, j ^L^ 


(27) 


where go is the long wave window channel . The clear column radiance 
L 

is computed from the equation 


I CLR,i,j (V) - [ Wi (V) - < 28 > 


which naturally follows from (27) once N.*. is defined. In this 

if 3 

method specific criteria are used to determine if FOVs are overcast, 
affected by broken clouds, or unaffected by clouds. No clear column 
computations are used for the overcast condition . Where cloud cor- 
rections are made as above (and in clear areas) an average value is 
obtained for the sub-grid area involved. 

(5) Chahine (1970) One-Layer Cloud Model . Assuming a single 
layer of clouds behaving as a black-body in equilibrium with the 
local ambient temperature, then cloud effects can be filtered from 
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the measured radiance if P^ (pressure at cloud top) and N (fraction 
of the FOV affected by clouds) can be specified. It can be shown 
that 


where 


I MEAS 


(V,P) 


W v ' p) - NG(v ' P c } 


G = I 


CLR 


- I 


CD 


(29) 


(30) 


The pressure at cloud-top height (P c ) may be obtained from one set 

of radiance measurements if the cloud— top temperature is known or 

from two sets of measurements made over adjacent areas with different 

cloud cover. Based on experimental evidence using simulated data, 

two different sets of radiance data may be used to determine P as 

c 

profiles obtained from the two sets of data should coincide from 
p c to the to P of the atmosphere (this concept is also basic to the 
Jastrow-Halem procedure that follows) . However, the equation 


N 


(I CLR " I MEAS )//G 


(31) 


a functional transform of the unknown temperature profile, should not 
then be used for some cloud-dependent frequency to obtain N from (31) 
If this result is used in the RTE for a partly cloudy atmosphere 
[Ref. (39)] by substituting N for A, then the corresponding residuals 


™ = IW V) - I,,) I / 'W V) <32) 

are small for all cloud-dependent frequencies for any temperature 

profile below clouds. Therefore, an extra parameter such as surface 

temperature is necessary to obtain N. Chahine concludes that any one 

of the combinations (P , N) , (P ,T ), (T ,N) , or (T , T ) will suffice 

c c o c co 

to allow determination of the temperature profile. Then assuming 

the combination (T ,T ) is known, (31) is used for a cloud-dependent 
c o 

frequency to calculate N and thus retrieve a new temperature profile 
from the RTE for a partly cloudy atmosphere. Using the retrieved 
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profile a new value of N is obtained from (31) and the process is 
repeated until convergence is obtained. A more detailed discussion 
of the RTE for a partly cloudy atmosphere will be presented in 
Section 3. 

(6) Jastrow and Halem (1973) . This procedure is a modification 
of Chahine's cloud model. Cloud height is first computed by deter- 
mining the initial near approach value (Ref. Chahine (1970) model) 
of the profiles obtained from sets of radiance observations obtained 
over partly cloudy areas but calculated assuming no clouds. Next, 
the radiance emitted over the cloudy portion of the ares is calcu- 
lated (computed from top of atmosphere to cloud-top height) . Then 
this computed radiance, I CD (v), is used in the equation 


I (v) = (1 - N) I (v) + NI (v) 
MEAS CLR CD 


(33) 


to calculate the equivalent clear column radiance, I (v) , for sev- 

CljR 

eral values of N between zero and one. By interpolation, the value 
of N is selected that produces the clear column radiance used to 
calculate the temperature profile whose surface temperature is 
closest to an observed surface temperature value . 

(7) Chahine (1974) . Assuming two FOVs with different amounts 
of cloud at the same height, clear column radiance values are computed 
from a guessed temperature profile for all frequencies used. Next, 
the clear column value for a "cloud dependent" window channel, w, is 
used in the equation 

n = [i to) - / n to) - WMl • < 34 > 

1 L ^ 


In this equation the subscripts denote the two FOVs and r\ is related 
to the fractional cloud cover by 


n = N x / (N 2 - n x ) 


( 35 ) 


Clear column radiance is then constructed from the expression 


20 



W v) “ + n[I ME AS 1 (V) - I MEAS 2 (V)1 > 0 • (36> 

A temperature profile is retrieved through use of this clear column 
radiance. All steps are then repeated to determine new values of 
computed clear column radiance/ I(o>), constructed clear column radi- 
ance (i.e. equivalent clear column radiance) , and the temperature 
profile. The iteration is continued until convergence of the clear 
column radiance values is attained. Chahine points out that as 
I^, T1D (v) is dependent on the temperature profile the problems of 
determining the cloud coefficient and the temperature profile are 
"...inseparable and should be carried out simultaneously." Other 
variations of the method presented above are given in the same ref- 
erence, including a single FOV, dual-frequency approach. The 
methods have been tested with simulated data and are based, as the 
1970 model, on Chahine' s relaxation scheme using the frequency set 
Vj to recover temperatures at solution points T(P _.) • and a cloud 
frequency (or frequencies for single FOV) to determine N. 

(8) Chahine (1975) . This model is unique in that temperature 
profiles are retrieved without calculating the clear column radiance. 
A single layer of black -body clouds is assumed, and an analytical 
transformation is derived to relate the temperature profile that 
would be derived from clear column radiance values directly to the 
apparent temperature profile. The apparent temperature profile is 
defined as the profile obtained from the measured radiance data with- 
out accounting for cloud effects. Illustrations are provided for 
simulated observations in the 15-ym band. For the single FOV, pro- 
files showed good agreement with the two FOV solutions obtained and 
the exact profile when the fractional cloud cover is less than five 
tenths, but not as good for increased amounts of cloud cover. 
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3. A NEW METHOD FOR RETRIEVING TEMPERATURE PROFILES FROM THE RADI- 
ATIVE TRANSFER EQUATION FOR A PARTLY CLOUDY ATMOSPHERE 

a. General 

In previous sections, methods for retrieving temperature pro- 
files from the RTE for a clear atmosphere have been discussed. In 
the method of Chahine (1968, 1970) a given frequency and a pressure 
level are paired to retrieve temperature solution points at the 
specified pressure levels. Smith (1970) presented an iterative 
equation in which there " ...is no limiting assumption made about 
the analytical form of the profile imposed by the number of radi- 
ance observations available." Duncan (1974a) demonstrated that 
Chahine' s relaxation formula could also be used with no limiting 

assumption imposed by the number of sounding frequencies utilized. 

* 

He then used Smith's method of calculating temperature values at 
each atmospheric level from a weighted average of the temperature 
values calculated for each frequency. 

In the discussion of cloud models used in profile retrievals it 
was mentioned that Chahine (1970) used the calculated cloud parame- 
ters with the RTE for a partly cloudy atmosphere to retrieve tempera- 
ture profiles. This method was based, as Chahine 's subsequent work, 
on the pairing of frequency with pressure level mentioned previously 
in order to achieve convergence over the widest possible spectrum of 
guessed temperature profiles. Now, a method will be presented which 
extends Duncan's procedures to an RTE for a partly cloudy atmosphere, 
and a new cloud model will be devised to accommodate observed cloud 
parameter (a) in this method. 

b. The method 

A solution to the problem of obtaining temperature profiles 
directly from cloud-attenuated radiance measurements would have to 
account for the heights, amounts, and opacities of the clouds that 
appear in the FOV. For a single layer of clouds assuming zero 
cloud reflectivity, Fritz ejt al. (1972) presented the equation 
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(37) 


I rn (v) = e (v)I (v) + [1 - e (v) 

CD C D O 


]i 


CLR 


(V) 


where is cloud emissivity and I^fv) is the radiance associated 
with black-body cloud conditions. By substitution of (37) into (33), 
the equation 


W (v) “ A V v > 


+ (1 - 


A >W V > 


(38) 


where A=Ne was obtained (assumes clouds are gray bodies) . 

We consequently write the RTE for a partly cloudy atmosphere as 

x. 


I(V) = A B[v,T(x p )]-t[v,y(Xp )] - f ° B [V , T (x) ] 

c c J 


9t [Vf y (x) ] 
9x 


dx 


(39) 


+ (1 - A) B[v,T(x ) ] * t[v,y (x 
o 


x 

°"-J 


° Brv,T(x) I- — 1 — d x 

dX 


where P^ is the cloud-top pressure . In an analogous manner to the 
clear case we may then define an expression 


B[v. ,T (n+l) (x) ] = X * (V i } B [v . ,T (n) (x) ] 
1 I*(v i ) 1 


(40) 


where I*(v.) is the measured cloud-contaminated radiance, and I*(\j.) 

i 

is the computed estimate of this unfiltered radiance. As the Planck 
function values are constant for a given frequency and level even 
though appearing under the integral sign, (40) will satisfy the 
RTE, and if the residuals 


R (n) (v.) = I*(v.) - I* (v. ) 

l li 


(41) 


are sufficiently small then the guessed temperature profile is a so- 
lution to the RTE. If the guessed temperature profile is not a solu- 
tion, T (n+1) pj ma y then be computed theoretically as in the 
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cloudless case from the weighted average 


,(n+l) 


N 


(x) =Y^ T (n+1) (V i# x)W(v ± ,x) / W(v ± ,x) 


N 


(42) 


i=l 


i=l 


where 


W(v ,x) = dtfv^x) if x ^ x^ or x p 


(43) 


W(v. , x ) = t(v. / x ) (44) 

1 O 1 o 

and 

W(v.,x p ) = At(v.,x p ) + (1 - A)dt(v ± ,x p ) # (45) 

1 c c 

Here the W values for each frequency are also the weighting functions 
used to compute radiance from the RTE for a partly cloudy atmosphere. 
However, Duncan (1974b) has demonstrated that for the atmospheric 
temperature range the Planck function can be approximated with suf- 
ficient accuracy by a Taylor's series expansion about a guessed 
temperature through the first derivative term. Although independent 
estimates of temperature for each frequency are not obtained, the 
approximations 

T (n+1) ( X ) = T (n) (x) + AB[T (n+1> <x)]/(.01T (n) (x) - 1.3) (46) 

where AB[T n+ ^ (x) ] is the weighted [weights are given in (43) - (45)] 
average of AB[v. ,T^ n+ ^ (x) ] , and 

B[V i ,T <n+1) (x)] = B[v i ,T (n+:L) (x) 1 + AB[T (n+1) (x)] (47) 

save significant machine computational time with no noticeable loss 
of accuracy. They have, therefore, been used throughout this 
research. 

In using (40) as a basis for calculating subsequent temperature 
values, the equation becomes less valid below clouds unless it is 
assumed that the ratio of the measured-to-calculated radiance values 
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is approximately equal to the ratio of calculated clear column radi- 
ance to the radiance that would have been measured in a clear situa- 
tion [Ref. (3)]. 

The assumption is required because the amount of measured radi- 
ance arising from below cloud level decreases as the amount of cloud 
cover increases. Therefore, as cloud cover increases, calculated 
temperatures below the cloud layer are increasingly based on 
radiances arising from cloud-top level. The assumption is also 
required because in calculating a weighted average [Ref. (42)], sig- 
nificant weight will frequently be given to channels that peak below 
cloud-top level. If a correct combination of cloud amount and height 
were used and the guessed and true profiles were equivalent in 
temperature at every level, then the assumption would be completely 
valid. As long as errors in the guessed profile do not get too big, 
errors in the retrieved profile should remain small even though 
cloud cover increases. If the assumption is generally valid, the 
degree of cloud contamination should not significantly affect tem- 
perature calculations below cloud level. Furthermore, if the ratio 
of measured to calculated radiance is greater (less) than one it 
can be seen from (40) that the temperature at each level of the 
guessed profile must increase (decrease) in order to provide a bet- 
ter estimate of measured radiance. Therefore, as cloud cover 
increases the importance of choosing an initial guessed profile that 
is in error in the same direction (with respect to sign) both above 
and below clouds assumes greater importance . As the shape of the 
guessed profile should not change significantly with successive 
iterations (Chow, 1975) , a smooth guessed profile beginning with a 
known surface temperature should* generally produce a temperature 
profile that is more accurate that the guessed profile even below 
clouds . 

c . The cloud model 

A temperature profile may therefore be retrieved through use of 
(39) and (40) if accurate values of effective cloud cover and cloud- 
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top height can be obtained. For multiple observed cloud layers, the 
equation 

IMv.) = £ A.(I b (v.,.) + (1 - A Total ) W v.> (4: 

j = l 

where k is the number of cloud levels, can be used. But these are 
black-body clouds (A=N) with the A values as would be observed from 
the top of the atmosphere looking down. It will be shown later that 
observed cloud cover is not an effective tool in retrieving tempera- 
ture profiles. However, an estimate of cloud height alone may be 
used as a first guess in calculating effective cloud cover for a 
single layer through use of the procedures described below: 

(1) Use the guessed temperature profile and only the window 
channel to calculate an estimate of the measured radiance for a 
clear situation from the RTE. Next, .subtract the calculated value 
of clear column radiance in the window channel from the measured 
value of radiance in the window channel. If the sign of the resul- 
tant value is positive, then the radiance calculated assuming no 
cloud has been found to be less than the measured radiance for a 
cloudy atmosphere. This result must, of course, be erroneous. The 
attempt to calculate effective cloud amount must therefore be aban- 
doned or a revised estimate of surface temperature employed in the 
calculations. However, if the sign of the resultant is negative, 
we may proceed to the next step. 

(2) Calculate radiance in the window channel for an overcast 
(A=1.0) situation at the estimated level and subtract this calcula- 
ted value of radiance from the measured radiance in the window 
channel. If the sign resulting from this subtraction is negative, 
then the radiance calculated from an overcast atmosphere is greater 
than the measured radiance. The estimated cloud-top temperature is 
therefore too warm. If the estimated cloud-top height is known to 
be accurate, then a revised cloud-top temperature may be used. How- 
ever, in this research reasonable estimates of the guessed profile 
are employed while estimated cloud-top height is assumed to be in 
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error. Therefore, overcast black-body cloud radiance is again calcu- 
lated and subtracted from the measured radiance at successively higher 
levels until a positive sign results from the subtraction performed. 
When a positive sign is obtained, whether it results from calcula- 
tions using the initially estimated height or from the trial and 
error method described above, we may proceed to calculate effective 
cloud cover at the chosen level . 

(3) In' the present research effective cloud cover was calculated 
through use of a one dimensional search. The search is performed at 
the chosen level as follows: 

(a) Calculate radiance for the window channel using a 
value of effective cloud cover that halves the possible choices for 
a value that will satisfy the RTE (i.e. first use 0.5) . 

(b) Subtract the value of radiance calculated in (a) from 
the measured radiance value. 

(c) If the sign of the results of the subtraction in (b) 
is negative, then the value we are seeking must lie between 0.5 and 

1 . 0 . 

(d) If the sign of the results of the subtraction in (b) 
is positive, then the value we are seeking lies between 0.0 and 0.5. 

(e) Calculate radiance for the window channel using a 
value of effective cloud cover that halves the possible choices for a 
value that will satisfy the RTE (i.e. either 0.25 or 0.75). Repeat 
the procedures discussed above (i.e. halve intervals) until the dif- 
ferences between measured and calculated radiance in the window 

channel approaches zero. 

2 

Smith (1976) has noted the one dimensional search described above 
is based on the equation 

I* (05) - A[I b (05)] + (1 - A) i clr (05) = 0 . (49) 

He therefore suggests that in order to save computer time the equation 
2 

Smith, W. L., 1976: Personal communication. 
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( 50 ) 


A 


i* (t o) - i cl ^M 

I b ( “> - W U) 


be used instead of the search. 

In any case, computed effective cloud amount is not unique except 
for a perfect guessed profile (guessed temperature = true value at all 
levels) when it is calculated at the true cloud level. However, re- 
sults that will be presented in the next section seem to indicate 
that if the estimated cloud-top height is not grossly in error, then 
choosing a cloud amount and height combination that gives the best 
estimate of measured radiance for the guessed profile will lead to a 
better retrieval than when one of the cloud parameters is correct 
but the other is in error. Also, the possibility of using the re- 
trieved profile to calculate an improved value of cloud cover and 
height and, hence, improve the retrieval will be investigated. 
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4 . PARAMETRIC STUDY 
a. Procedures 

Prior to using the method described in the previous section with 
real data it was tested with simulated radiance data (radiances cal- 
culated through vertical integration of the RTE from a known tempera- 
ture profile and known cloud parameters) to determine the effect of 
errors in observed or calculated parameters on the retrieved tempera- 
ture profile. The general procedure was to cause an error in one or 
more of the parameters needed to retrieve the temperature profile 
while other parameters retained their true values (i.e. the values 
used to calculate synthetic radiance data) . The method outlined in 
Section 3b is particularly well suited to such an approach as the 
effects of varying cloud parameters are immediately apparent in the 
retrieved profile since cloud parameters are not filtered out before 
commencing the retrieval. 

Simulated radiance data were calculated for the channels of the 
NOAA-2 satellite through use of (39) or (48) and the temperature pro- 
files designated "true" in the figures to follow. The fractional 
cloud layers and their heights used to calculate the simulated radi- 
ances are noted at the top of each figure and also adjacent to the 
pressure value corresponding to cloud height at the side of each 
diagram. 

In the temperature profile retrieval, the guessed profile used 
is a version of the 1962 U.S. Standard Atmosphere with 10 °C added to 
the temperature at each level above the surface through 221 mb. The 
guessed profile is then approximately five degrees cooler than the 
true profile for most of the troposphere at all levels except the 
surface. Cloud parameters that are read in are designated "input 
cloud" at the top of the figure while those that are calculated are 
designated "calc cld." Profiles are shown up to 250 mb. 

As in the Duncan model, numerical integration for the model was 
accomplished through use of the trapezoidal rule using centered dif- 
ferences for a maximum of 100 pressure levels from 1000 mb to 0.01 mb 
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that are equally spaced in pressure to the two-sevenths power. Where 
surface pressure falls between two levels the higher level is used 
as the surface for computational purposes (e.g. surface temperature) 
except that the centered difference at the bottom of the atmosphere 
was applied between the surface pressure and the pressure level above 
the next highest pressure level. 

b. Results 

In Figs. 2 and 3 the effects of not accounting for clouds in the 
presence of a poor surface temperature guess are illustrated. As 
anticipated from the work of Chahine (1970) , in each case the apparent 
temperature profile begins to coincide with the true profiles in the 
vicinity of the top of the cloud layers . Figures 4 and 5 (same 
guessed profile and cloud parameters as Figs„ 2 and 3, respectively) 
are illustrations of the improvement made in the retrieved profiles 
through use of the procedures described in Section 3c to calculate 
cloud cover when the cloud-top height is known. Surprisingly, im- 
provement was noted even in profiles where very large amounts of 
cloud cover were present (Fig. 4) . This result is probably due to 
the validity of the assumption discussed in the previous section for 
the NOAA-2 weighting functions and the guessed profile used. Tests 
of profiles calculated for several values of fractional cloud cover 
and height reveal that errors of ±1 level (=1000 ft in the tropo- 
sphere) have a relatively insignificant effect on the retrieved pro- 
file. This magnitude of error is probably representative of those 
that would be made in handling actual data for many low cloud situa- 
tions. The effect of making a larger error in a high cloud situation 
will be illustrated later in this section. 

Profiles were retrieved for several cloud parameter values from 
guessed profiles exhibiting sharp inversions. Results of a typical 
retrieval are illustrated in Fig. 6. There is no doubt that the 
calculated profile takes on the shape of the guessed profile. In the 
absence of prior knowledge, the necessity of using a smooth guessed 
profile in connection with this model becomes obvious. 
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Pressure (mb) 


i 


EFFECTIVE CLOUD =0.80 
HEIGHT = 1200m AGL 
ASSUMED CLEAR IN RETRIEVAL 



Fig. 2. Temperature sounding retrieved from Duncan's Method for 
a clear atmosphere using simulated radiance measurements 
for broken (0.80) low cloud conditions. There is a +5°C 
error in the guessed surface temperature. 
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EFFECTIVE CLOUD COVER =0.18 

HEIGHT = 8500m AGL 

ASSUMED CLOUDLESS IN RETRIEVAL 



-20 -10 0 10 20 30 

Temperature °C 

Fig. 3. Temperature sounding retrieved from Duncan's Method for 


a clear atmosphere using simulated radiance measurements 
for scattered (0.18) high cloud conditions. There is a 
+5°C error in the guessed surface temperature. 



Pressure (mb) 


EFFECTIVE CLOUD =0.80 
HEIGHT = 1200m AGL 
CALC CLD = .52 



Fig. 4. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere. Simulated radiance measurements were 
prepared for broken (0.80) low cloud conditions. There 
is a +5°C error in the guessed surface temperature. 
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EFFECTIVE CLOUD = 0.18 
HEIGHT = 8500m AGL 
CALC CLD = .23 
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5. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere. Simulated radiance measurements were 
prepared for scattered (0.18) high cloud conditions. 

There is a +5°C error in the guessed surface temperature. 
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Pressure (mb) 


EFFECTIVE CLOUD =0.60 

HEIGHT = 2200m AGL 

CALC CLOUD = 0.69 at 2200m 



Fig. 6. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing a guessed 
profile exhibiting sharp temperature inversions. The 
simulated radiance measurements were prepared for broken 
(0.60) middle level cloud conditions using a +5°C error 
in the guessed surface temperature. 
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Testing indicates that observed cloud parameters could be used 
to retrieve profiles even for the two-layer case (Fig. 7) if their 
values could be specified exactly. However, small errors in these 
parameters lead to unacceptable retrieval errors (Fig. 8) . 

An attempt was made to retrieve profiles from simulated data de- 
rived from two-layer cloud cases through use of the one-layer cloud 
model. Here a single cloud layer that returns radiance values that 
approximate those of the true cloud parameters is sought. Examples 
of results are shown in Figs . 9 and 10 . In these examples the esti- 
mated single cloud-layer height was chosen as the highest cloud layer. 

The accuracy of retrieved profiles deteriorates when the height 
of a significantly lower layer is chosen as the estimated height for 
a single layer. As calculated radiance for the low layer must be 
larger than the measured radiance, a higher level is sought to calcu- 
late the fractional cloud amount [Ref. Paragraph 3c- (2)]. Since this 
process does not give unique values of fractional cloud amount/height 
it should be noted that the greatest amount of calculated fractional 
cloud cover will occur at the lowest possible height. In Fig. 11 
an example of the results of choosing a significantly low estimate 
of average cloud-top height is shown. 

An example of a profile retrieved through use of a guess that is 
significantly too high is given in Fig. 12 where the estimated height 
is 1900 m higher than the true highest cloud tops. It should be 
noted that estimating cloud-top heights too high will give radiance 
values that are lower than comparable measured values . Since calcu- 
lated clear column radiance values will be higher than the measured 
radiance, the sign change required by the model is accomplished and it 
is possible to calculate a radiance that is approximately equal to 
the true radiance at all levels above the true cloud tops (assuming 
one-layer cloud cover) . For more than one layer an estimated level 
somewhat higher than the highest average cloud-top layer might be 
appropriate (Ref. Fig. 10 where calculated cloud amount is greater 
than the true amount at the highest level) . 

Finally, numerous attempts were made to improve the accuracy of 
retrieved profiles by calculating an initial cloud amount from the 
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EFFECTIVE CLOUD = 0.40+0.20 
HEIGHTS = LOW( 3200m AGL) , HIGH(9500m) 
CLD WAS INPUT AND NOT CALCULATED 



Fig. 7. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing in the re- 
trieval the same values of the cloud parameters as used 
to calculate the measured (simulated) radiance values. 

The simulated radiance measurements were prepared for 
two layers of clouds with the fractional amount and height 
values shown above. 
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EFFECTIVE CLOUD = 0.50 + 0.30 
HEIGHTS = LOW (3200m AGL) , HIGH(9500m) 

INPUT CLOUD =0.40 at 3200m + 0.20 at 9500m 



Fig. 8. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing in the re- 
trieval the same cloud-top heights as used to calculate 
the measured (simulated) radiance values but fractional 
cloud amounts that are each 0.1 less than the values used 
to calculate the simulated radiances . The simulated radi- 
ance measurements were prepared for two layers of clouds 
using the effective cloud amount and height values 
shown above . 
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EFFECTIVE CLOUD COVER = 0.40 + 0.20 
HEIGHTS = LOW(2200m AGL) , HIGH(8500m) 
CALC CLD =0.27 at 8500m AGL 



-20 -10 0 10 20 30 


Temperature °C 

Fig. 9. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing in the re- 
trieval the highest cloud-top height used to calculate 
simulated radiance measurements and calculating a one 
level fractional cloud amount (0.27) 'at that level. The 
simulated radiance measurements were prepared for two 
layers of clouds using the effective cloud amount and 
height values shown above . 
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EFFECTIVE CLOUD COVER - 0.50+0.50 
HEIGHTS = LOW (22 00m AGL) , HIGH(8500m) 
CALC CLD = 0.59 at 8500m AGL 



-20 -10 0 10 20 30 

Temperature °C 


Fig. 10. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing in the re- 
trieval the highest cloud-top height used to calculate 
simulated radiance measurements and calculating a one 
level fractional cloud amount (0.59) at that level. The 
simulated radiance measurements were prepared for two 
layers of clouds using the effective cloud amount and 
height values shown above . 
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EFFECTIVE CLOUD = 0.50+0.50 
HEIGHTS = LOW ( 2200m AGL) , HIGH(8500m) 
CALC CLD = 0.95 at 5000m 



-20 


-10 


0 


10 


20 


30 


-20 -10 0 10 20 30 


Temperature °C 

Fig. 11. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing in the re- 
trieval a cloud height that is significantly lower than 
the highest cloud-top height used to calculate simulated 
radiance measurements and calculating a one level frac- 
tional cloud amount at the significantly low level. The 
simulated radiance measurements were prepared for two 
layers of clouds using the effective cloud amount and 
height values shown above . 
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EFFECTIVE CLOUD = 0.40+0.20 
HEIGHT = LOW {2200m AGL) r HIGH(8900m) 
ESTIMATED HEIGHT - 10,800m AGL 
CALC CLD - 0.18 at 10,800m AGL 



-20 -10 0 10 20 30 


Temperature °C 


Fig. 12. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere accomplished by employing in the re- 
trieval a cloud-top height that is significantly higher 
than the highest cloud-top height used to calculate sim- 
ulated radiance measurements and calculating a one level 
fractional cloud amount at the significantly high level. 
The simulated radiance measurements were prepared for two 
layers of clouds using the effective cloud amount and 
height values shown above. 
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guessed profile as previously, but then using the retrieved profile 
to calculate a new estimate of cloud amount and continuing this 
iteration until measured and calculated radiance values converged 
for all channels. These attempts generally resulted in slightly 
improved profiles. Figure 13 is an example of a profile obtained 
in this manner. However, this procedure is unsatisfactory when 
applied directly to real as opposed to simulated data. 
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EFFECTIVE CLOUD = 0.80 
HEIGHT = 2200m AGL 
ESTIMATED HEIGHT = 2500m AGL 

CALC CLD = 0.25 0.55 (CONV IN 4 ITERATIONS) 



-20 -10 0 10 20 30 

Temperature °C 


Fig. 13. Temperature sounding retrieved from the RTE for a partly 
cloudy atmosphere. It was accomplished by calculating 
a value of fractional cloud amount at an estimated cloud- 
top height, employing the calculated value to retrieve 
a temperature profile and then calculating revised values 
of cloud amount and temperature until calculated radiance 
values and the simulated measurements converge. The sim- 
ulated measurements were prepared for two layers of 
clouds using the effective cloud amount and height values 
shown above. 
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5. APPLICATION TO ATMOSPHERIC VARIABILITY EXPERIMENT (AVE III) DATA 

a . Procedures 

In Section 3 a method was developed to retrieve temperature pro- 
files directly from cloud- contaminated radiance data. In Section 4 
the method was examined through use of simulated data. However, these 
data are easier to handle and less noisy than real data. Furthermore, 
transmittance errors do not affect simulated data. The feasibility 
of using the model for a real case was therefore investigated. The 
general procedure followed was to use NOAA-4 satellite radiance data 
measured over the area of the AVE III experiment to retrieve tempera- 
ture profiles, and to compare the retrieved profiles with the excellent 
radiosonde data obtained during the AVE III experiment for the sta- 
tion nearest the center of the applicable radiance spot (resolution 
element) . Highest cloud-top heights and surface temperature were 
estimated from the synoptic observations available for the time clos- 
est to satellite passage over the area. No attempt was made to refine 
the cloud-top or. temperature data from any other data source or to 
correct the shelter temperature in any way to give a closer approxi- 
mation of the true surface value for the area of the radiance spot. 

b. Data 

A description of the Vertical Temperature Profile Radiometer 
(VTPR) used on NOAA series satellites and the accuracy and format of 
retrieved data is given by McMillin el: al . (1973) . The VTPR scans 

from left to right in 23 discrete steps per scan line. In the 0.5-sec 
interval allowed for each spot, radiance measurements are obtained 
in six channels of the CO^ band, a window channel, and a channel in 
the water vapor absorption band (not used in this study) . Incremental 
steps of about 2.7° are used, giving 30.3° areal coverage from the 
nadir direction for each scan line. When viewing in the nadir di- 
rection the projection of each scan spot on Earth's surface is approx- 
imately a square 55 km on a side, with spot size increasing somewhat 
with viewing angle. After completing a scan line (12.5 sec of which 
11.5 is used to make measurements) the instrument takes an additional 


45 



second to return to its original position. Scan spots are contiguous 
along and across the satellite track. The orientation of the eight 
lines of the radiance data used in this experiment with reference to 
Earth's surface and AVE III observing stations in the area is shown 
in Fig. 14. Scanning was initiated at approximately 0233 GMT, 

6 February 1975. 

AVE III radiosonde data for 6 February 1975 at approximately 
0000 GMT were used for comparison with the temperature profiles re- 
trieved from satellite radiance data. The method of collecting and 
processing the AVE III data was discussed by Fuelberg and Turner 
(1975) . Sounding data were obtained for 51 stations (Fig. 15) , 11 
in the area of interest of this study, for every pressure contact 
and interpolated to give values at 25-mb intervals from Earth's sur- 
face to 25 mb. RMS errors for computed temperature values are esti- 
mated to be 1 1°C . As temperatures are required for approximately 
100 intervals from 0.01 mb to the surface in the present study, the 
necessary values were obtained through linear interpolation of the 
AVE III temperatures. In like manner dew point values were obtained 
to 258 mb, the highest level for which these data are used in the 
program to compute weighting functions . Surface temperatures and 
dew points were obtained from surface synoptic observations or in- 
terpolation from 0000 GMT radiosonde surface values . Temperature 
and dew point values for Stephenville, Texas, were not available above 
375 mb, therefore, measurements in this region of the atmosphere for 
Shreveport; Louisiana, were substituted for the missing Stephenville 
data. The substitution was necessary for use in computer programs 
to compare the AVE III data with the satellite-derived profiles, but 
the comparisons for Stephenville above 375 mb are, of course, invalid. 

The surface weather map for 6 February 1975 at 0000 GMT (after 
Fuelberg and Turner) is shown in Fig. 16. The data in Table 1 were 
used to estimate a cloud-top height for each station. These heights 
were then used as input to the cloud model previously described. Al- 
though there is no known way to retrieve the temperature profile under 
an overcast layer of black-body clouds from the radiation originating 
from the surface under these clouds, retrievals were attempted for 
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Fig. 15. 


Rawinsonde stations participating in the AVE III 
experiment (Fuelberg and Turner, 1975) . 
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Table 1. Estimation of average cloud-top heights 


Station 

SFC P (mb) 
(From AVE III) 

Cloud and Visibility 
Observations 
02Z 03Z 

Estimated P (mb) 
at Cloud Top 

Centerville, Ala 

996.6 

30 010 

30 010 

870 

Jackson, Miss. 

1003.0 

280- CD 20+ 

280-0 20+ 

329 

Shreveport, La. 

1008.5 

Not Available 

870 

S tephenvi lie, Tx 

971. -2 

Ml 6 015 

M18V ©15 

870 

Del Rio, Tx. 

977.7 

250-0 20+ 

250-020+ 

361 

Midland, Tx. 

913.3 

250-© 20+ 

250 020 

344 

Nashville, Tenn. 

989.4 

M2 40 65 010 

240 M35© 10 

699 

Little Rock, Ark 

1007.1 

M2 20 12 

M22 012 

902 

Monet te , Mo . 

966.1 

Ml 10 10 

Mil © 10 

810 

Amarillo, Tx. 

889.3 

60M35© 10 

1O0M35012 

699 

Marshall Space 

Flight Center, Ala. 991.2 

700 250 015 

2500 15 

377 
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all stations including those where overcast conditions were observed. 
It was felt that at the overcast stations, as long as a relatively 
smooth guessed profile was used, the assumption discussed in Section 
3c, i.e. 


t*(V.) I(v ± ) 

__ ^ I (v ^y 


(51) 


would be valid and an improvement of the guessed temperature values 
below cloud level could be obtained. It is also possible that al- 
though a station reports an overcast condition, it may not be a true 
black -body overcast. The overcast layer may be thin, and may not 
completely attenuate radiance arising from below the cloud layer. 
Because of its low emissivity, a reported overcast layer of thin 
clouds may yield no greater value of effective cloud cover (A=Ne) than 
a broken or even scattered layer of thick clouds . 


c. Results 

In order to retrieve an accurate temperature profile through use 
of the method outlined in this study it is necessary to employ a first 
guess profile that approximates the true temperature values . Both 
smooth climatological profiles and numerical forecasts have been used 
with the "minimum information" method as guessed profiles. The 
question of the proper guessed profile to choose for a particular 
retrieval is necessarily dependent on the knowledge of the true pro- 
file that is possessed prior to attempting the retrieval. As can be 
seen from Fig. 6, the shape of an inversion in the guessed profile 
will be picked up in the retrieved profile and at approximately the 
same atmospheric pressure . The coincidence of shape thus attained is 
independent of the validity of the inversion. Thus to attain the 
most accurate retrieval it would seem best not to include fine detail 
in the guessed profile. However, if a retrieval is desired for a 
location at the same latitude and in the same air mass as a station 
for which a recent radiosonde run is available, it would seem logical 
to use at least some of the known radiosonde temperatures (with mod- 
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ification in the lower layers and smoothing between radiosonde tem- 
peratures as necessary) as the guessed profile in the retrieval. Even 
so, for situations comparable to the 6 February 1975 situation, where 
an area is under the direct influence of a strong frontal system, 
using the radiosonde run of one station as the guessed profile of 
a nearby station might lead to fictitious features in the retrieved 
profile. This would also be true if a day old profile for a station 
were used as the guessed profile. Furthermore, it was desired that 
guessed profiles for all stations be coincident in temperature above 
the lower portion of the atmosphere so that changes in similar pro- 
files might be noted for the various stations. 

In investigating the 6 February situation several different 
guessed profiles were tried. Guessed profiles were computed as fol- 
lows: (1) Temperatures were averaged at each level above the surface 

to 25 mb for the 51 stations of the AVE III experiment. The profile 

3 

of temperature above 25 mb was provided by Dr. L. Duncan (1975) and 
is a standard profile for the White Sands area. (2) Temperatures 
were averaged as above for the 51 AVE III stations. Next, a constant 
lapse rate between levels was computed from the surface to level 90 
(699 mb) and a new constant lapse rate for each 10 levels thereafter 
to level 50 (97 mb) was computed. These values were substituted 
for the previously computed averages below level 50. (3) The same 

method as in (2) above was employed except averages were computed 
for only the eleven stations for which radiance data were obtained. 

It was hoped that this procedure would lead to a smooth profile that 
incorporated the general shafce of the average profile. (4) Tempera- 
tures were computed as in (3) , but the dew point profile was computed 
differently. For the three profiles discussed above the 0000 GMT 
AVE III dew points for each station were employed. In this case dew 
points were averaged and smoothed in precisely the same manner as the 
temperature profiles in (2) , (3) , and (4) . 

For each of the guessed profiles discussed above retrievals were 
performed for each of the eleven AVE III stations . The root mean 

3 

Duncan, L. D., 1975: Personal communication. 
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square errors from the surface to cloud-top level resulting from com- 
parison of the true and retrieved profiles were calculated. Results 
are shown in Table 2. Chahine (1970) has noted that the effects of 
clouds on the retrieved profile are almost entirely confined to the 
region below cloud-top level. This effect may also be noted in Figs. 

2 and 3. Comparison of RMS errors in Table 2 for the same stations 
but different guessed profiles should therefore provide a relative 
estimate of the usefulness of the various guessed profiles tested. 
Results shown in Table 2 indicate that in guessed profile (1) , a 
simple average of the true temperature profiles at all AVE III sta- 
tions/ false detail in the guess causes serious errors in the re- 
trieved profiles. Comparing results of profiles (1) and (2) , it can 
be seen that smoothing out the false detail will invariably result in 
an improved retrieval. Use of the eleven-station average points in 
profile (3) does not give significantly better results than retrievals 
obtained using guessed profile (2) . Neither does use of an average, 
but smoothed, dew point profile [profile (4) ] significantly improve 
the below cloud-top level retrieved profiles . All changes in the 
guessed profile have little effect on the large RMS errors calculated 
for Stephenville and Midland. The errors are probably the result 
of a gross error in estimated cloud-top height caused by the presence 
of clouds that are much higher than estimated. Improvement of re- 
trievals at these stations will be attempted later in this section. 
From Table 2, guessed profile (4) appears to give the best results 
and is used for subsequent retrievals unless otherwise indicated. 

In the retrievals discussed above a convergence interval of 0.1 
2 -1 

mw/m sr cm was used for the convergence of computed to measured 

radiance values for each channel . After the computation of the cloud 

parameters was achieved a total, of fifty iterations was allowed to 

achieve convergence in the interval prescribed. Both the interval 

and the maximum number of iterations had been used in computations 
4 

by Duncan (1975) . However, for the 6 February 1975 cases convergence 
was in no case achieved within 50 iterations for all channels. 


4 


Ibid. 
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Table 2. RMS errors (°c) below cloud-top level for various guessed 
profiles . 


Station 

(1) 

Retrieved 

vs 

True 

(2) 

Retrieved 

vs 

True 

(3) 

Retrieved 

vs 

True 

(4) 

Retrieved 

vs 

True 

Centerville 

6.8 

3.8 

2.9 

2.3 

Jackson 

3.2 

2.3 

* 2.0 

2.0 

Shreveport 

0.9 

1.4 

1.3 

1.6 

Stephenville 

* 

4.8 

5.0 

4.7 

Del Rio 

4.5 

2.9 

2.7 

2.1 

Midland 

* 

4.3 

4.2 

4.1 

Nashville 

2.6 

1.5 

1.6 

1.8 

Little Rock 

1.4 

2.0 

1.5 

1.7 

Monet te 

7.5 

3.4 

3.5 

3.3 

Amarillo 

3.7 

1.7 

1.9 

2.1 

MSFC 

* 

3.0 

2.9 

3.0 


♦Profile not computed. 
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Increasing the convergence interval beyond the error tolerance allowed 
in computing cloud cover is impractical as convergence does not occur 
at the same rate in each channel. As pointed out by Chahine (1970), 
the relaxation method upon which the present research is based is 
"...a discrete numerical process in which the concept of formal con- 
vergence plays hardly any role..." and the rate of convergence is 
judged as the rate at which the residuals for each channel reach their 
"asymptotic" values. For the profiles investigated, decreasing the 
number of iterations below fifty revealed that differences between 
measured and computed radiance reached different, but nearly constant 
values, for the various channels in a fairly rapid manner. A review 
of results (RMS errors) obtained by comparing AVE III profiles with 
profiles retrieved using various numbers of iterations (Tables 3 and 
4) gives the impression that the degree of accuracy can rarely be 
improved after approximately ten iterations . This is in agreement 
with the empirical results found in tests of Chahine" s method (Conrath 
and Revah, 1972) . Comparison of RMS errors in the guessed profile 
with RMS errors in retrieved profiles (Table 4) reveals that acqurate 
retrievals for ten iterations were generally obtained to tropopause 
level. With the exception of Marshall Space Flight Center (MSFC) , 
errors greater than those that are present in the guessed profile 
generally occur only when inversions are present in the guessed or 
true profiles. At MSFC the guessed profile fits the true atmosphere 
so well that iteration does not appear to lead to improvement of the 
guessed profile. However, for ten iterations the maximum RMS error 
in the results is only 1.0 degree greater than the guessed value 
(Table 4) . 

It was assumed that results for some of the overcast and high 
cloud cases could have been degraded by a gross error in the estimate 
of the height of the highest cloud layer in the field of view of the 
radiometer. Therefore, all profiles were recomputed using an estimated 
cloud-top height of 299 mb. Significant improvement was noted at 
Stephenville and Midland, and the computed cloud parameters were used 
therefore in subsequent research . Results for ten and one iterations 
are shown in Table 4 . 
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Table 3. RMS errors (°C) for retrieved profiles. 


Station 

SFC 

RMS 

to CLD 

error (°C) for 
SFC to 726mb 

iterations shown 
699 to 489mb 

469 

to 314mb 


75 

50 

25 

75 

50 

25 

75 

50 

25 

75 

50 

25 

Centerville 

2.6 

2.3 

2.2 

2.7 

2.1 

1.9 

2.3 

1.7 

1.5 

4.9 

4.1 

3.2 

Jackson 

2.6 

2.0 

1.5 

2.7 

2.3 

2.1 

1.9 

1.3 

0.9 

3.0 

2.2 

1.2 

Shreveport 

2.0 

1.6 

1.5 

2.2 

1.6 

1.4 

3.0 

2.9 

2.9 

4.2 

3.5 

2.6 

Stephenville 

4.6 

4.7 

4.7 

3.4 

3.6 

3.9 

8.3 

7.9 

7.0 

8.0 

7.3 

5.9 

Del Rio 

2.8" 

2.1 

1.6 

2.3 

1.5 

0.8 

2.4 

1.7 

1.5 

3.7 

2.9 

2.2 

Midland 

4.1 

4.1 

4.2 

2.8 

3.0 

3.3 

4.2 

4.5 

4.8 

4.8 

4.3 

4.9 

Nashville 

1.9 

1.8 

1.8 

1.0 

1.0 

1.0 

2.1 

2.0 

2.3 

4.1 

3.3 

1.9 

Little Rock 

1.7 

1.7 

1.6 

1.6 

1.6 

1.6 

3.8 

3.3 

2.3 

7.0 

6.2 

4.9 

Monette 

3.8 

3.3 

2.9 

3.7 

3.0 

2.3 

3.8 

3 .1 

2.7 

1.8 

0.7 

2.2 

Amarillo 

2.1 

2.1 

2.0 

1.9 

1.9 

1.8 

1.9 

1.7 

1.8 

4.3 

3.7 

2.7 

MSFC 

3.5 

3.0 

2.4 

3.9 

3.4 

2.8 

3.6 

3.1 

2.6 , 

2.2 

1.5 

0.8 


*True profile for Shreveport used above 375 mb. 
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**True profile for Shreveport used above 375 mb 



In Figs. 17 through 27, the guessed, retrieved (using 10 itera- 
tions) , and AVE III profiles are shown to 300 mb. It should be noted 
that in all cases the errors in the surface temperatures were probably 
greater than computed because the observed shelter temperatures were 
used as both the true and first guess temperatures. However, errors 
for the first few levels above the surface would normally be less 
than shown as the AVE III temperatures at these levels were not mod- 
ified for the time lag between the times the radiosonde and satellite 
measurements were taken. 

Isotherms were drawn for the nearest level to standard from 850 
to 300 mb for the AVE III and retrieved temperature data (Figs. 28- 
37) . Surface frontal positions where shown are from the analysis by 
Fuelberg and Turner (1975) . 

Comparison of the analyses shows the most serious error in grad- 
ient resulted from the spurious warm ridge of Fig. 31. This feature 
was caused by the very poor 699-mb temperature at Monette. Thus the 
primary reason for this occurrence may be traced to the fictitious 
inversion used in the guessed profile for Monette where the 699-mb 
surface is at the apex of the inversion. The importance of not intro- 
ducing fictitious features in the guessed profile (in this case the 
inversion was created by using average data across a frontal zone) is 
again emphasized. Analysis of the retrieved isotherms in the lower 
levels appears consistent with surface frontal positions. In eval- 
uating the retrieved analyses it should be recalled that at all levels 
above 839 mb guessed temperatures were the same for each station at 
a particular level. Improvement was thus demonstrated, even in a 
generally cloudy situation. Greater improvement would be anticipated 
when account is taken of anticipated profile features for a given air 
mass or station. 

An exact comparison of the results achieved for the cases shown 
in Table 4 and the isotherm maps with the results of previous investi- 
gators is impossible. This is not only because of the different 
methods of retrieval and different guessed profiles employed, but also 
the variance in three-dimensional space investigated, the averaging 
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Fig. 17. Centerville, Ala. 6 Feb 1975 retrieved temperature pro- 
file compared with the guessed profile and AVE III 
radiosonde data. 
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Fig. 18. Jackson, Miss. 6 Feb 1975 retrieved temperature profile 
compared with the guessed profile and AVE III radio- 
sonde data. 
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Fig. 19. Shreveport, La. 6 Feb 1975 retrieved temperature pro- 
file compared with the guessed profile and AVE III 
radiosonde data. 


Guessed Profile 
Retrieved (0233Z) 
AVE III (OOOOZ) 


— a — 



61 




Pressure (mb) 
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Fig. 20. Stephenville, Tx. 6 Feb 1975 retrieved temperature pro- 
file compared with the guessed profile and AVE III 
radiosonde data. 
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Fig. 21. Del Rio, Tx. 6 Feb 1975 retrieved temperature profile 
compared with the guessed profile and AVE III radio- 
sonde data. 
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Fig. 22. Midland, Tx. 6 Feb 1975 retrieved temperature profile 
compared with the guessed profile and AVE III radio- 
sonde data. 
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Fig. 23. 
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Nashville , Tenn. 6 Feb 1975 retrieved temperature pro- 
file compared with the guessed profile and AVE III 
radiosonde data. 
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Fig. 24. Little Rock, Ark. 6 Feb 1975 retrieved temperature 

profile compared with the guessed profile and AVE III 
radiosonde data. 
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Fig. 25. Monette, Mo. 6 Feb 1975 retrieved temperature profile 
compared with the guessed profile and AVE III radio- 
sonde data. 
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Fig. 26. Amarillo, Tx. 6 Feb 1975 retrieved temperature profile 
compared with the guessed profile and AVE III radio- 
sonde data. 
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(03Z) 250015 
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Fig. 27. Marshall Space Flight Center, Ala. retrieved temperature 
profile compared with the guessed profile and AVE III 
radiosonde data. 
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Fig. 29. Retrieved 839 mb isotherms and surface frontal positions. 
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of clear and cloudy results and the frequent reporting of results ob- 
tained exclusively from simulated data, and the comparison of results 
achieved with analyzed charts as well as radiosonde data from nearby 
stations. There are no reports of single FOV profiles achieved for 
exclusively cloud-contaminated cases in recent literature. The re- 
sults achieved in this study are considered significant in that 
improvement over a guessed profile was achieved for a single FOV. 

The improvement shown from surface to cloud top is especially signi- 
ficant because it is in this region that profiles retrieved through 
use of previously-mentioned single FOV cloud models tend to deteriorate 
as cloud cover increases . 

The results achieved for overcast cases require explanation as 
these cases cannot normally be handled through use of previously pub- 
lished cloud models. 

For tne model presented in the present study a low estimate of 
the top of a single layer overcast requires a search at successively 
higher levels until computed radiance from the estimated overcast top 
to the top of the atmosphere is less than measured radiance. The 
lowest level in the atmosphere where this result is possible must 
correspond to the level of greatest cloud cover possible of all the 
possible cloud amount and height combinations which will yield the 
correct radiances . With an accurate estimate of the top of the over- 
cast/ retrieval of a reasonably accurate profile from this point to 
the top of the atmosphere should be possible regardless of the re- 
trieval method. Where thick clouds are present above the overcast/ as 
was apparently the case at Stephenville, a poor estimate of the known 
overcast height will give a poor retrieval. 

It appears that when cloud conditions were correctly specified/ 
retrieved temperature profiles below the overcast layer also showed 
soma improvement compared to the guessed profile even though the con- 
tribution of the atmosphere below cloud level to the measured radiance 
values was assumed to be nil. It is believed that the improvement 
occurred because the Planck functions upon which temperatures at all 
levels are based were determined through use of the ratio of measured 
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to computed radiance [Ref. (40)], regardless of the contributions of 
specific levels to the total radiance values., For any one iteration, 
computed temperatures at all levels in the troposphere should increase 
(decrease) if all measured radiance values arc greater (less) than 
computed radiance values for channels whose weighting functions peak 
in the troposphere. This result would be anticipated if guessed 
profile temperatures were less (more) at all levels than the true 
atmospheric temperature values. The initial guessed profiles used 
in the present study usually approximated this situation in that they 
were uniformly greater or uniformly less for a given location than 
true values at corresponding levels. As retrieved profiles tend to 
retain the shape of the guessed profile (i.e. adjustments to the 
guessed profile occur nearly uniformly with height) , the relationship 
of true to calculated temperatures below an overcast tended to show 
improvement . 

d. Applicability to mesometeoro logical research 

It is believed that the procedures employed in the present re- 
search may have applicability to mesometeorological research. For any 
time at which a resonable guessed profile can be forecast and radi- 
ance data is available, temperature profiles can be retrieved and 
used to fill the data gaps between the synoptic hours. Also, excellent 
spatial resolution is achieved through use of a single FOV method, and 
it appears feasible to study patterns of temperature change over 
relatively small areas through use of retrieved temperature profiles. 

It appears that the best procedure to follow in retrieving profiles 
for mesometeorological research would be to use a guessed profile ob- 
tained through use of a known shelter temperature, a few other widely 
separated tropospheric temperature values obtained by averaging over 
an air mass or a section of an air mass, and a constant lapse rate 
between the chosen temperature values. Where fine detail of the true 
atmospheric profile is thought to be present, it would appear best in 
most cases to use a relatively smooth profile as a guess and then 
add the suspected detail to the retrieved profile. 
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e. Possible sources of error 

No attempt was made to provide an accurate estimate of the tem- 
perature profile above the troposphere. To accomplish this would 
require that an estimate of the tropopause height be included in the 
guessed profile for each station. The problem is essentially the 
same as encountered when attempting to account for other inversions 
in the true profile. If structure such as the tropopause is intro- 
duced in the first guess, it must be present or serious errors are 
caused (Wolski, 1975) . However, it is also to be expected that fail- 
ure to introduce the tropopause, and consequently its correct contri- 
bution to the calculated radiance values, would lead to some error 
in the retrieved profiles at all levels. 

Another possible source of error is the procedure for calculating 
the weighted average used in temperature computations. 

An attempt was made to modify the retrieval computer program so 
that temperature at a given atmospheric level in the troposphere was 
computed through use of (40) but only at the frequency that provided 
maximum input to the measured radiance at that level (e.g. channel 
4 was used between 500 and 300 mb [Ref. Fig. lb]) . Retrieved profiles 
were uniformly less accurate than for comparable retrievals using 
previously-discussed procedures and displayed significant computational 
instability . 

The question of whether a special weighting function [Ref. (45) ] 
is required in computing the weighted average used for temperature 
calculations was also investigated. Equation (43) was therefore used 
instead of (45) in temperature calculations (but not to calculate 
radiance values) and the resulting retrievals were compared with those 
obtained through use of (45) . Significant differences were only 
noted for the overcast cases. At Monette, comparison with the AVE III 
profile revealed that errors in the retrieved temperatures using (43) 
were over 2K more than most comparable retrieved values using (45) . 
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6 . 


CONCLUSIONS 


A method has been presented to retrieve single FOV tropospheric 
temperature profiles directly from cloud- contaminated radiance data 
through use of auxiliary data such as observed shelter temperatures 
and estimated cloud-top height.. The iterative technique utilized was 
an extension of the work of Chahine (1970) as modified by Smith (1970) 
and Duncan (1974a) . A model was formulated to calculate cloud para- 
meters for use with the RTE through use of a one-dimensional search 
[or (51)3 at an estimated cloud-top level where it has been shown to 
be possible to calculate an effective cloud amount that will satisfy 
the RTE and provide an approximation of measured radiance for the 
guessed profile. 

The method was evaluated through use of simulated data and for 
a coincident data sample from the AVE III experiment and NOAA-4 satel- 
lite for an area dominated by an active cold front and covered by 
considerable cloudiness at various levels. 

The major conclusions derived from the present research are: 

(1) A single FOV method of retrieving temperature profiles from 
cloud-contaminated radiance data that improves the accuracy of guessed 
profiles has been developed. Through use of a single FOV method many 
temperature profiles can be retrieved for the same area in which a 
single average temperature profile can be retrieved through use of 

a multiple FOV technique. It is significant that in the method 
presented improvement in the guessed profile was noted under the 
cloud layer where retrievals using other single FOV techniques tend 
to deteriorate. The method requires estimates of surface temperature 
and average cloud- top height that are not grossly in error. 

(2) It is possible to make an accurate estimate of the average 
tops of a thick overcast layer through use of the cloud model developed 
as an integral part of the single FOV retrieval method discussed above 
when there are no thick clouds present above the overcast layer and 
the guessed temperature profile is relatively accurate . 

(3) For most overcast situations it should be possible to achieve 
accurate retrievals at least down to cloud-top level. 
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(4) Observed cloud parameters are not obtained with sufficient 
precision to use directly in the RTE. Significant errors in retrieved 
profiles resulted when this procedure was tested. 

Though not conclusions, the following items should be noted: 

(1) Through judicious choice of a guessed profile, it appears 
possible to improve guessed profiles independent of the cloud amount 
present. Profiles retrieved during the parametric study from guessed 
profiles that were uniformly colder or warmer at all levels than the 
true values exhibited this characteristic. 

(2) Improvement in the guessed profile through utilization of 
the procedures discussed in this study should occur whenever a reason- 
able estimate of the true lapse rate is forecast. However, the abso- 
lute accuracy of the retrieved profile is also a function of the 

c* priori knowledge of the state of the atmosphere possessed by the 
researcher. 

(3) It appears that suspected detail should not be included in 
the guessed profile, but might profitably be added to the retrieved 
profile . 

(4) Use of the techniques described above to provide useful 
data for mesometeorological research appears feasible for any time 
radiance data is available and a reasonable guessed profile can be 
forecast. A guessed profile utilizing a constant lapse rate between 
average air mass temperature values known with some accuracy {e.g. 
obtained from NMC analyses) should normally lead to an accurate 
retrieval. Also, through use of the many temperature profiles that 
can be retrieved over a relatively small horizontal area with a single 
FOV method, it should be possible to determine an accurate pattern 

of temperature change. 

(5) The largest errors in retrieved profiles should be antici- 
pated in the vicinity of moving frontal disturbances. Knowledge of 
the observed shelter temperature reduces this source of error. 

(6) Due to the fact that the observed shelter temperature was 
used in the guessed profiles, fictitious features were introduced when 
the upper-level guess was based on data averaged across a frontal zone. 
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(7) No significant improvement in retrieved profiles from 
simulated data was obtained by refining estimates of cloud cover 
through use of successively retrieved profiles . For AVE III data 
deterioration in retrieved profiles was noted. 
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7. 


RECOMMENDATIONS FOR FURTHER RESEARCH 


The following specific suggestions are presented: 

(1) Further testing of the method outlined in the present re- 
search for various cloud conditions and guessed profiles would provide 
useful information . 

(2) The extent of degradation of the retrieved temperature pro- 
files caused by errors in tropopause height should be determined. 

(3) The method should be tested with data from the Nimbus series 
of satellites . 

(4) The extent of improvement of retrieved profiles that could 
be obtained through use of microwave data should be determined. 

(5) Retrieved profiles could be compared with profiles obtained 
using various cloud models and/or retrieval techniques . 

(6) Smith et_ al^. (1974) presented a method for determining 
effective cloud height which may be used with auxiliary data for a 
single FOV. It would be interesting to compare results obtained 
through use of these procedures with those obtained using the proce- 
dures discussed above for the same data set. This suggestion was made 
by Smith (1976) 5 . 

(7) Using available radiosonde data at surrounding points, tem- 
perature profiles should be estimated. The extent that these profiles 
can be improved through use of the method outlined in the present 
research should be examined. 

(8) A temperature profile estimated from surrounding radiosonde 
data should be used to prepare a guessed profile for a specific point. 
If surface temperatures and cloud observations are available, then 
temperature profiles may be retrieved for a relatively small area; and 
the pattern of temperature change , the thermal wind and temperature 
gradients may be examined. 

(9) In the cloud model outlined in the present research an 
estimated height level was rejected if calculated window channel radi- 
ances did not lie between radiances calculated for zero and overcast 

5 Smith, W. L., 1976: Personal communication. 
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(1.0) cloud cover for the same level. These limits could be refined 
for many cases to provide better estimates of cloud height. For 
example, if 0.3 of cloud is observed, limits of 0.1 and 0.5 might be 
tested. 

(10) As mentioned previously, no significant improvement in 
retrieved profiles -was obtained by refining estimates of cloud cover 
through use of successively retrieved profiles. However, both for 
the simulated and AVE III data, this procedure was tested while allow- 
ing 50 iterations on a guessed profile once a cloud amount had been 
estimated. It may be that through use of 10 or less iterations com- 
putational instability could be avoided and better results achieved. 
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